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The existence and properties of the gravitational phase transition in this model are in- 
vestigated numerically. It is shown that, by cooling a nondegenerate gas of weakly inter- 
acting massive fermions below some critical temperature, a condensed phase emerges, 
consisting of quasidegenerate fermion stars. For fermion masses of 10 to 25 keV, these 
fermion stars may very well provide an alternative explanation for the supermassive 
compact dark objects that are observed at galactic centers. 
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Abstract 

The Thomas-Fermi model at finite temperature is extended to describe a system of 
self-gravitating weakly interacting massive fermions in a general-relativistic framework. 



1 Introduction 

The ground state of a condensed cloud of fermionic matter, interacting only gravitationally 
and having a mass M below the Oppenheimer-Volkoff (OV) limit [[|], is a cold fermion 
star in which the degeneracy pressure balances the gravitational attraction of the fermions. 
Degenerate stars of fermions in the mass range between 10 and 25 keV are particularly 
interesting 0, as they could explain, without resorting to the black- hole hypothesis, at least 
some of the features observed around the supermassive compact dark objects with masses 
in the range of M — 10 6 ' 5 to 10 9 ' 5 solar masses, that are reported to exist at the centers of a 
number of galaxies [§, f|, [|, ||, [S], j8j, including our own || |1(J, and quasistellar objects (QSO) 
111, |12 , 13| , |i~4f . Indeed, there is little difference between a supermassive black hole and a 



fermion star of the same mass near the OV limit, a few Schwarzschild radii away from the 
object jE^, . 



The purpose of this paper is to study, in the framework of a general-relativistic Thomas- 
Fermi model, the formation of such a star that could have taken place in the early universe 
shortly after the nonrelativistic fermionic matter began to dominate the radiation. This 
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system was previously studied in the Newtonian approximation [O, [TR [19], 12(1 ETL |22|, and 



it was shown that the canonical and grand-canonical ensembles for such a system have a 
nontrivial thermodynamical limit |I7|, Under certain conditions these systems will 

undergo a phase transition that is accompanied by a gravitational collapse [^TJ |2^|. The phase 
transition occurs uniquely in the case of the attractive gravitational interaction of neutral 
fermions. As the phase transition does not happen for particles obeying Bose-Einstein or 
Boltzmann statistics, this phenomenon is quite distinct from the usual gravitational clus- 
tering of collisionless dark-matter particles. Gravitational condensation will also take place 
if the fermions have an additional short-range weak interaction, as neutrinos, neutralinos, 
gravitinos, and other weakly interacting massive particles do. 

Effects of general relativity cannot be neglected when the total mass of the system is close 
to the OV limit 0. There are three main features that distinguish the general- relativistic 
Thomas-Fermi model from the Newtonian one: i) the equation of state is relativistic, ii) 
the temperature and chemical potential are metric-dependent local quantities, and in) the 
gravitational potential satisfies Einstein's field equations instead of Poisson's equation. 

This paper is organized as follows: In Section we briefly discuss the nonrelativistic 
Thomas-Fermi model at finite temperature. In Section ^ this model is extended within 
a general-relativistic framework. In Section [| we discuss the solution at zero and finite 
temperature and, in particular, the conditions under which the first-order gravitational phase 
transition occurs. Conclusions are drawn in Section |5| and, finally, in Appendix |A], we prove 
a theorem on the extremal properties of the free energy. 



2 Thomas- Fermi model in Newtonian gravity 

Consider a system of N fermions of mass m interacting only gravitationally, confined in a 
spherical cavity of radius R, in equilibrium at a finite temperature T. For large N, we can 
assume that the fermions move in a spherically symmetric mean-field potential V(r) which 
satisfies Poisson's equation 

-4n(rV) = 4irGm 2 n, (2.1) 
r dr 2 

G being the gravitational constant. The number density of the fermions (including an- 
tifermions) n can be expressed in terms of the Fermi-Dirac distribution (in units H = c = 
k = l) 

( g 2 V(r) yi"- 1 



n{r) = g 



d 3 q 



1 + eXP l v 2^T + T 



(2.2) 



(27T) 

Here g denotes the combined spin-degeneracy factors of the neutral fermions and antifermions, 
i.e., g is 2 or 4 for Majorana or Dirac fermions, respectively. For each solution V(r) of 
Eq. ( |2.1|) , the chemical potential \i is adjusted so that the constraint 



J d 3 rn(r) = N (2.3) 
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is satisfied. It may be shown that a particular spherically symmetric configuration n(r) will 
satisfy Eqs. (|2.1|)-(|2.3|) if and only if it extremizes the free energy functional defined as |18 



F[n] = fi[n]N - - J d 3 rn(r)V[n] 
f d 3 rd s q ( 

- T9 lw ( 1+exp 



2mT 




(2.4) 



where V[n] and p[n] are implicit functionals of n(r) through Eqs. ( p.l|) - (p.3|) . For a physical 
solution, we have to require that the free energy is minimal, i.e., 
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> 0. 



(2.5) 



The set of self-consistency equations (|2.1| )- (|2.3| ), together with (|2.5| ), comprises the nonrela- 
tivistic gravitational Thomas-Fermi model. 

It may be easily shown that the following scaling property holds: If the potential en- 
ergy V(r) is a solution to the self-consistency equations (|2.1| )-( f2~3l) , then the rescaled V = 
A^V(Ar), with A > 0, is also a solution with the rescaled fermion number N = A 3 N, radius 
R = A~ l R, and temperature T = A A T. This property, which will be referred to as nonrela- 
tivistic scaling, implies the existence of a thermodynamic limit of N~ 7 ^ 3 F, with A^ 1//3 i? and 
N~ 4 I 3 T approaching constant values for iV — > oo. In this limit, the Thomas- Fermi equation 
becomes exact fTTR 



3 Thomas- Fermi model in general relativity 

Consider a self-gravitating gas consisting of iV fermions of mass m in equilibrium within 
a sphere of radius R. We denote by p, p, n, and a the pressure, energy density, particle 
number density, and entropy density of the gas, respectively. The metric generated by the 
mass distribution is static, spherically symmetric, and asymptotically flat, i.e., 

ds 2 = i 2 dt 2 - (1 - 2M/r)- 1 dr 2 - r 2 (d0 2 + sinfld0 2 ). (3.1) 

Einstein's field equations are then given by 

d£ = M + 4nr 3 p 

dr 5 r(r - 2M) ' 1 ' ' 



dM , 



dr 



4vrr>, (3.3) 



with the boundary conditions 



«fl) = (i-^y /2 ; M{0)=0 (34) 
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The equation of state may be represented in a parametric form [23| 



f d 3 Q E 

P = 9 J pi^W' (3 - 6) 

p=aT Swf ln(1 + e " £/f+A/f ) > ( 3 - 7 ) 

where (7 denotes the spin degeneracy factor and = \/m 2 ~+ g 2 . The quantities T and // 
are the local temperature and chemical potential, respectively. As discussed in Appendix [A], 
thermodynamic and hydrostatic equilibrium in the presence of gravity implies 

The constants T and p are the temperature and chemical potential at infinity. Although 
the matter is absent at r = 00, the temperature at infinity has a physical meaning: T 
is the "red shifted" temperature of the black-body radiation of a gravitating object in 
equilibrium at finite temperature measured at infinity. In other words, our gravitating object 
is in equilibrium with a heat bath which could be thought of as a black-body radiation in 
the empty space surrounding the object. As a consequence of Eq. ( |3.8| ), different gravitating 
configurations of the same size "in contact" with the same heat bath may have different 
surface temperatures. Therefore, the relevant equilibrium parameter is the temperature of 
the heat bath, T, and not the surface temperature of the gravitating object. Particle number 
conservation yields the constraint 

R 

dr 4vrr 2 (l - 2M /r)~ 1/2 n{r) = N. (3.9) 







Given the temperature at infinity T, the set of self-consistency equations (|3T^ ) -( fO|) defines 



the general-relativistic Thomas-Fermi equation. One additional important requirement is 
that a solution to the equations (f3 . 2|) - (pT9|) should minimize the free energy. Based on the 
considerations given in Appendix [A|, the free energy may be written in the form 

F = M + fxN- J dr4vrr 2 £(l - 2M/r)- 1/2 (p + p), (3.10) 

with M = M.(R). The theorem proved in Appendix A guarantees that solutions to Eqs. 
( |3.2| )-( |3~9D extremize the free energy F, i.e., the free-energy functional assumes either a 
maximum or a minimum. We only have to find out which of the solutions are maxima and 
discard them as unphysical. 

Next we show that, in the Newtonian limit, we recover the usual Thomas- Fermi model 
as defined in Section 0. Introducing the nonrelativistic chemical potential unr = \i — m 
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and the approximations £ = e v3 ~l + y?, -E~m + q 2 /2m and Ai/r -C 1 , we arrive at the 



Thomas- Fermi self-consistency equations [21, 22 in the form 



P f d?q ( q 2 m fi NR ,\ 1 , . 

n = ^ = 9 J j2^[ +eM ^ + T v --r ) ) ' (3 - n) 

dip M dM A _ 2 



dr r 2 ' dr 



Anr 2 p, (3.12) 



mN 

<pm = ~-\ -M(0) = 0, (3.13) 

R 

drA-Kr 2 n(r) = N, (3.14) 

o 



which are equivalent to the set of equations fl2~T| ) - ([2.3|) . The free energy (|3.10|) in the New- 
tonian limit yields 

1 r R r R 

F = mN + PnrN / dr Anr 2 rup — / dr4irr 2 p, (3.15) 

2 Jo Jo 

with 

-^/IH 1 ^-^-^^)- (316) 



which, up to a constant, is equal the nonrelativistic Thomas- Fermi free energy 

A straightforward thermodynamic limit iV — > oo, as discussed by Hertel, Thirring, and 



Narnhofer |18|, pLOfl ■» is not directly applicable in the general-relativistic case. First, in con- 
trast to the Newtonian case, there exists a limiting configuration with maximal M and N 
(the OV limit) at zero temperature, and, as we shall shortly demonstrate, also at finite 
temperature. Second, the scaling property of the relativistic Thomas-Fermi equation, which 
will be referred to as relativistic scaling, is quite distinct from nonrelativistic scaling. This 
scaling property may be formulated as follows: If the configuration {£( r ) ? (r)} is a solu- 
tion to the self-consistency equations (|3 . 2|) - ( pT9|) , then the configuration {£ = ^(A^r), M. = 
AJ\A.(A~ l r); A > 0} is also a solution with the rescaled fermion number iV = A 3 ^ 2 N, radius 
R = AR, asymptotic temperature T = A~ l l 2 T, and fermion mass m = A~ x l 2 m. The free 
energy is then rescaled as F = AF. Hence, there exists a thermodynamic limit of N~ 2 ^ 3 F, 
with N~ 2 / 3 R, N l / 3 T, and N l l 3 m approaching constant values when iV — > oo. 



4 Numerical integration 



In the following we use the units in which G = 1. We choose appropriate length, mass and 
fermion number scales a, b, and c respectively, such that 



or, restoring h, c, and G, we have 



^hMpx ^ oior [2 /l5keV\ 
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1.3185 x 10 10 W- km, (4.2) 

V g cm 2 V 9 \ m J 

[2 Ml m [2 Zl5keV\ 2 , . 

b=J ^ = 0.8929 x 10 10 W- M 4.3 

y g m z y g \ m J 

c =. pJk = 5.5942 X 10", P- f«' (4.4) 

where Mpi = J%cj~G denotes the Planck and M the solar mass. 

We are looking for a solution of the Thomas- Fermi problem as a function of temperature. 
For numerical convenience, let us introduce a new parameter 

« = f (4-5) 

and the substitution 

£= ii($ + i)-V2. (4. 6 ) 
m 

Using this and (|4.1|) , Eqs. (|3.5|) - (|3.7|) may be written in the form 

1 r 00 y 2 

n = - dy £ , (4.7) 

- 1 + eMUiv 2 + !)/($ + !) " l)a] 
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P = -» / / — > (4.^ 

" - 1 + exp[(J(y2 + !)/($ + 1) - l)a] 



JO 

1 r 00 w 4 
p = / dy z (4 g) 

3n 2 Jo v ^TT(1 + exp[(^/(y 2 + l)/($ + 1) - l)a] ' 

In this way, both the fermion mass and the chemical potential are eliminated from the 
equation of state. 

The field equations (CO) and (|3.3|) now read 



d<& M. + Anr 3 v 

^ = -2($ + l Y y*™* , (4.10) 
dr v ' r{r- 2M) V ' 

dM , 2 



dr 



Anr 2 p. (4.11) 



To these two equations we add 



dM 
dr 



47rr 2 (l - 2M /r)- 1/2 n (4.12) 
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imposing the particle-number constraint as a condition at the boundary 



Af(R) = N. (4.13) 
Equations ( f4.10|) -( [4.12|) should be integrated using the boundary conditions at the origin 

$(0) = $ > -1 ; A^(0) = 0; ; W(0) = 0. (4.14) 
The parameter $o, which is uniquely related to the central density and pressure, will even- 



tually be fixed by the requirement fl4.13|) . For r > R, the function $ yields the usual 
empty-space Schwarzschild solution 

^ r) = 4( 1 _?5£)- 1 _ 1 , (4,5) 



with 

M = M(R) = / dr47rr 2 p(r). (4.16) 
Jo 

We now show that a solution to the general-relativistic Thomas-Fermi equation exists 
provided the number of fermions is smaller than a certain number N max that depends on a 
and R. From ( |4.8| ) and ( [4.9| ) it follows that, for any a > 0, the equation of state p{p) is an 
infinitely smooth function and dp/ dp > for p > 0. Then, as shown by Rendall and Schmidt 
|25fl , there exists for any value of the central density p a unique static, spherically symmetric 



solution of the field equations with p — > as r tends to infinity. In that limit Ai(r) — > oo 
and Af(r) — > oo, as may be easily seen by analyzing the r — > oo limit of Eqs. fl4.10|) and 
( |4.11| ). However, the enclosed mass M and the number of fermions N within a given radius 
R will be finite. We can then cut off the matter from R to infinity and join the interior 
solution onto the empty space Schwarzschild exterior solution by making use of equation 
( |4.15| ). This equation together with ( |4.5| ) fixes the chemical potential and the temperature 
at infinity. Furthermore, it may be shown that our equation of state obeys asymptotically 
at high densities a 7-law, i.e., p = const n 1 and p = (7 — l)p, with 7 = 4/3. In this case, 
as is well known [£6], there exists a limiting configuration {?/' 00 (r), A4(r)oo} such that M 



and iV approach nonzero values ¥«, and iVoo, respectively, as the central density p tends 
to infinity. Thus, the quantity iV is a continuous function of po 011 the interval < po < 00, 
with iV = for p = 0, and iV = as p — > 00. The range of N depends on a and R and 
its upper bound may be denoted by N max (R, a). Thus, for given a, R and iV < A^" max (_R, a) 
the set of self-consistency equations (fOD-( f06l) has at least one solution. 

As is evident from the equation of state (|3.5|) - (|3.7|) , if we do not fix the boundary and 
do not constrain the particle number N, the pressure (and the density) will never vanish 
(except perhaps at r = 00), unless T = 0. Thus, since we fix the boundary at R and cut 
off the matter from R to infinity, the pressure (and the density) will have a discontinuity. 
This characteristic of the non-relativistic Thomas- Fermi model in atomic physics pTfl , and 



Newtonian gravity |18|, remains also in general relativity. However, the density and 

the pressure decrease rapidly with r, so if R is chosen sufficiently large, the pressure and 
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the density at the boundary will be extremely small. Furthermore, the region r > R is 
never empty in reality, so that a positive pressure at the boundary is more realistic than a 
vanishing pressure. 

The numerical procedure is now straightforward. For a fixed, arbitrarily chosen a, we first 
integrate equations (|4.10|) and (|4.11|) numerically on the interval (0, R) and find solutions for 
various initial $ - Simultaneously integrating ( |4.12j ), we obtain J\f(R) as a function of $ - 
The specific value of $o is then determined such that Af(R) = N. The chemical potential 
fi corresponding to this particular solution is given by ( f4.15|) . If we now eliminate /x using 
( |Q>1) , we finally get the parametric dependence on temperature through a. 

Let us first discuss a degenerate fermion gas (T = 0) as a reference point that can also 
be compared with the well-known results by Oppenheimer and Volkoff [0]. In this case, the 
Fermi distribution in fl4.7p -( fOD becomes a step function that yields an elementary integral 
with the upper limit yp = v$ related to the Fermi momentum 



Qf 



m 



The equation of state can be expressed in terms of elementary functions of $ 

1 



3tt 2 



P 



8tt 2 
1 

247T 2 



(2$ + 1 W$($ + 1) - ArshV$ 



(2$ - 3W + 1) + 3Arsh\/$ 



(4.17) 



(4.18) 

(4.19) 

(4.20) 
At this 



The radius of the star is naturally defined as the point where the density vanishes, 
point, owing to ( 4.18 ), $ = 0. Therefore, we integrate equations ( |4.7|) -( |4~9D starting from 
r = up to the point R where $(-R) =0. As a result, the quantities M, N, and R are 
obtained as functions of the parameter $ , which is related to the central particle-number 
density through ( |4.18| ). In Fig. [1|, we plot the mass of the star as a function of the radius 
R. The maximum of the curve corresponds to the OV limit p]. The limiting values are 
Roy = 3.357, M Q v = 0.38426, and N Q y = 0.39853, in units of a, b, and b/m respectively. 
The curve left from the maximum represents unstable configurations that curl up around 
the point corresponding to the infinite central density limit. 

We now turn to the study of nonzero temperature. The quantities T, N, and, R are free 
parameters in our model and their range and choice are dictated by physics. The temperature 
T is restricted only to positive values. The number of fermions iV is restricted by the 
Oppenheimer- Volkoff limit. The radius R is theoretically unlimited; practically, it should 
not exceed the order of interstellar distances. It is known that a classical, semi degenerate, 
isothermal configuration has no natural boundary in contrast to the degenerate case of zero 
temperature, where for given iV (up to the Oppenheimer- Volkoff limit) the radius R is 
naturally fixed by the condition of vanishing pressure and density. At nonzero temperature, 
if we, e.g., fix only N and T and let R — > oo, our gas will occupy the entire space, and hence 
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Figure 1: Mass versus radius for fermion stars at zero temperature in the general- relativistic 
framework (solid line) compared with the corresponding Newtonian approximation (dotted 
line). Dashed line is the black hole limit M = R/2. 



p and p will vanish everywhere. If we do not restrict iV and integrate the equations on the 
interval (0, oo), M and N will diverge at oo. In such a case one has to introduce a cutoff. In 



the isothermal model of a similar kind by Chau, Lake, and Stone ||30||, a cutoff was chosen at 



the radius R, where the energy density was about six orders smaller than the central value. 
Our choice is based on the following considerations: As in the Newtonian Thomas-Fermi 
model we expect that, for given number of fermions N < Nov, there exists a unique 
configuration that is a solution to the self-consistency equations Q3.5| ) to ( |3.9| ) and which 
becomes a degenerate Fermi gas at T = 0. For such a configuration, an effective radius 
R e s > -Rov m ay be defined so that $(i2 e ff) = 0. Although the density does not vanish at 
this point, most of the mass will be contained inside the sphere of radius _R e fr- If we choose 
the boundary at R 3> -Rov, the total mass will be dominated by the density distribution 
within R e g, and it will be almost independent of the choice of R. Thus, in the following we 
will work with R = 100 ~ SOifov- 

In Fig. |2|, the fermion number N is plotted as a function of initial $o for several values of 
the parameter a. In contrast to the T = case, all curves with finite a have a peak around 
$ = 0. The second peak corresponds to the OV limit. From this figure we can deduce that, 
for a given N, there is a range of a's for which the Thomas-Fermi equation may have more 
than one solution. This is a clear indication for the existence of an instability even below 
the OV limit and, as a consequence, we expect that a first-order phase transition occurs. 

Fixing N = 0.38, which is slightly below the OV limit, we can now plot the temperature 
as a function of a in Fig. S. Using this figure as a parametric function for temperature, 
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Figure 2: Fermion number N versus central potential $ f° r a — f — 300 (full line) 500 
(dashed line) and 150 (dot-dashed line). T=0 is represented by dotted line. 




Figure 3: Temperature T (in units of m) versus 1/a for N = 0.38 and R = 100. 
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Figure 4: Total mass (in units of b) versus temperature for N = 0.38 and R = 100. 

the mass, free energy, and entropy are shown as functions of temperature in Figs. |], [5], and 
respectively. In the temperature interval T = (0.0015 — 0.007)m there are three distinct 
solutions of which only two are physical, namely, those for which the free energy assumes a 
minimum. The solution that can be continuously extended to any temperature above the 
mentioned interval is referred to as "gas", whereas the solution that continues to exist at 
low temperatures, and eventually becomes a degenerate Fermi gas at T = 0, will be called 
"condensate" . In Fig. |2] the gas is represented on each curve by the part left from the first 
maximum, while the part from the first minimum up to the second maximum represents the 
condensate. By noting that $ is negative for the gas and positive for the condensate, we 
may define an order parameter as 

<f = $ + |$ |, (4-21) 

which is strictly positive in the condensed phase (ordered phase) and equal to zero in the 
gaseous phase (disordered phase). 

The phase transition takes place at the temperature T c , where the free energy of the gas 
and that of the condensate become equal. The dashed curves in Figs. [|, and ^| represent the 
physically unstable solution. In our example, the transition temperature is T c = 0.0043951m, 
as indicated in the plots by the dotted line. The latent heat per fermion released during the 
phase transition is given by the mass difference at the point of discontinuity 

= 0.0438m. (4.22) 

N 

So far, we have studied, as an example, an object with number of fermions N just below 
the OV limit iVov- Any object with N < iVov will undergo a gravitational phase transition 
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Figure 5: Free energy per fermion F/N (in units of m) versus temperature T for N = 0.38 
and R = 100. 



7 . 5 



2.5 




0.002 



0.006 0.008 



Figure 6: Entropy per fermion S/N versus temperature T for iV = 0.38 and R = 100. 
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Figure 7: Critical temperature as a function of the fermion number in the Newtonian (dashed 
line) and general- relativistic Thomas- Fermi model (solid line). 

at a critical temperature which depends on the mass, of course. With decreasing N, the 
cavity radius R must be appropriately increased, since the effective radius of the condensate 
increases following approximately the zero-temperature mass-radius relation. As iV becomes 
smaller, the system approaches the nonrelativistic scaling regime discussed in Sec. |2|, and 
for N <C iVov the critical temperature will decrease according to 

T c = const iV 4/3 , (4.23) 

if the cavity radius R is simultaneously rescaled as iV _1//3 . In Fig. [7] we compare the critical 
temperature calculated in both Newtonian and general-relativistic Thomas-Fermi models, 
as a function of N. The nonrelativistic scaling law turns out to be very accurate for iV < 
0.2 Ar ov . 

It is important to check that the critical temperature is not very sensitive to variations of 
the cavity radius R for the following two reasons: First, in our model R is arbitrary except for 
the requirement that it should be much larger then the effective radius which for N = 0.38 
is of the order R eS ~ _R Q v — 3.357. Second, if the critical temperature rapidly decreases 
with R, the adiabatic cooling of the gas through the universe expansion may not necessarily 
lead to the point of the phase transition. Fig. |8| shows that the critical temperature indeed 
decreases very slowly, roughly by a factor of two if R increases from 30 to 300. This is much 
weaker than the adiabatic cooling of a nonrelativistic gas which goes approximately as 1/R 2 . 
Thus we conclude that the gravitational phase transition will necessarily take place in the 
course of the universe expansion. 
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Figure 8: Critical temperature T c as a function of the cavity radius R for N = 0.38. 



5 Conclusions 



In this work, we have extended the Thomas-Fermi model to a general-relativistic framework. 
This model was then applied to a system of self-gravitating fermions. We have investigated 
numerically the circumstances under which this system undergoes a gravitational phase 
transition at nonzero temperature. This phase transition is quite distinct from the more 
extensively investigated strong-interaction driven phase transition that might occur in neu- 
~fl. The main underlying physics here is the competition between the partial 
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degeneracy pressure due to the Fermi-Dirac statistics and the attractive force due to the 
gravitational interaction. It is obvious that the application of this model to astrophysical 
systems will work very well if the non-gravitational interactions between the individual parti- 
cles can be neglected. This is certainly the case, e.g., for weakly interacting quasidegenerate 

I], [32], but perhaps it could be 



heavy neutrino neutralino, or gravitino matter [15, 22. |3 



valid even for collisionless stellar systems |33|, |34[] 



Finally, let us briefly comment on a similar model by Chau, Lake and Stone |30[ which 
was considered earlier in the context of a possible galactic massive neutrino halo. Their 
model differs from our approach in essentially two aspects: First, the equation of state is 
not consistent with the condition of thermodynamical and chemical equilibrium, i.e. with 
Eq. ( |3.8[ ) and second, the particle number constraint ( |3.9| ) is not imposed. Thus, in contrast 
to the Thomas-Fermi model discussed here, the Chau et al. model does not describe a 
canonical system in equilibrium. 
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A Free Energy 

Consider a canonical ensemble that describes a nonrotating fluid in equilibrium at nonzero 
temperature. We denote by w M , p, p, n, and a the velocity, pressure, energy density, particle 
number density, and entropy density of the fluid. A canonical ensemble is subject to the 
constraint that the total number of particles 

nu tM dE„ = N (A.l) 

should be fixed. The spacelike hypersurface £ that contains the fluid is orthogonal to the 
time-translation Killing vector field fc^ which is related to the velocity of the fluid by 

fcA» = £ u f. £=(A^) 1/2 - (A.2) 

The energy-momentum tensor is defined as 

T^ u = (p + p)*V** - pg^. (A. 3) 

The metric generated by the mass distribution in equilibrium is static, spherically symmetric, 
and asymptotically flat, i.e., 

ds 2 = i 2 dt 2 - X 2 dr 2 - r 2 (d9 2 + sin Odcf) . (A.4) 

The metric coefficients £ and A may be determined from Einstein's field equations. However, 
we shall for the moment assume that £ is an arbitrary function of r, which asymptotically 
approaches unity, and we parameterize A in terms of the mass M. within the radius r 

^(x-^lMp (A,) 

where r 

M(r) = [ dr' Anr' 2 p(r') . (A.6) 
Jo 

The temperature T and chemical potential p are, in general, metric-dependent local 
quantities. The local law of energy- momentum conservation, T^ u . u = 0, for a perfect fluid in 
a static gravitational field, yields the equation of hydrostatic equilibrium 



d„p = -(p + p)r 1 d v £- (A.7) 



15 



This equation, together with the thermodynamic identity (Gibbs-Duhem relation) 

4 = -4 + -4> ( A - 8 ) 

1 n 1 n 1 
and the condition that the heat flow and diffusion vanish 

| = const., (A.9) 

implies 

n = T- fe = n, (A.10) 

where T and \x are constants equal to the temperature and chemical potential at infinity The 
temperature T may be chosen arbitrarily as the temperature of the heatbath. The quantity 
H in a canonical ensemble is an implicit functional of £ owing to the constraint (|A.1|) . The 
first equation in (|A.1Q ) is the well known- Tolman condition for thermal equilibrium in a 
gravitational field [[36)1 . 

Following Gibbons and Hawking [57|, we postulate that the free energy of the canonical 
ensemble is 

F = M — [ TcrWdE^, (A.ll) 



where M is the total mass as measured from infinity. The entropy density of a relativistic 
fluid may be expressed as 

a= i(p + p-fxn), (A.12) 

T and /2 being the local temperature and chemical potential as defined in ( A.lOj ). Based on 
Eq. ( [A.10Q , the free energy may be written in the form analogous to ordinary thermodynamics 

F = M-TS, (A.13) 

with M = A4(R), and the total entropy S defined as 

S= [ R dr4irr 2 \l(p + p)-^N } (A.14) 
Jo 1 1 

where we have used spherical symmetry to replace the proper volume integral as 

/ u M dS„= [ R drAnr 2 X. (A. 15) 

Jt, Jo 

The following theorem relates the extrema of the free energy with the solutions of Ein- 
stein's field equation: 
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Theorem: Among all momentarily static, spherically symmetric configurations {^{r), M(r)} 
which, for a given temperature T at infinity, contain a specified number of particles 

i-R 

Aixr 2 dr X(r)n(r) = N (A.16) 



within a spherical volume of a given radius R, those and only those configurations that 
extremize the quantity F defined by ( |A.13|) will satisfy Einstein's field equation 



d£ M + 47rr 3 p 
dr r(r - 2M) 



with the boundary condition 



2M\ 1/2 



(A.17) 



(A.18) 



Proof: By making use of the identity ( |A.8| ) and the fact that S(p/T) = 5(p/T) and that N 
is fixed by the constraint (|A.16|) , from Eqs. ( |A.13| ) and (|A.14|) we find 



r R T r R T 

SF = 8M - dr4nr 2 7 =;(p + p)5X - drAnr 2 X^5p. 

Jo T Jo T 



(A.19) 



The variations 5X and 5p can be expressed in terms of the variation 5M(r) and its derivative 

d5M 



yielding 



dr 



47ir 2 5p, 
dX 



rR t dX r R 

5F = 8M- drA-nr 2 ^(p + p)——5M- dr 
Jo T oJvt Jo 



T5M 
T dr 

By partial integration of the last term, and replacing T/T by £, we find 

^r 2 i{p + p)^--Uxi) 



6F=[1- \(R)£(R)] SM - [ R dr 

Jo 



dM dr 



5M 



(A.20) 
(A.21) 

(A.22) 



where 5j\4(r) is an arbitrary variation on the interval [0, R], except for the constraint 
5.M(0) = 0. Therefore, 5F will vanish if and only if 



47rr 2 C(p + p) 



dX d 



dM dr 



(AO = 



(A.23) 



and 



1 - \(R)£(R) = 0. (A.24) 

Using Eqs. (|A~5| ) and (|A~6|) , we can write Eq. (|A.23|) in the form (|A.17|) , and Eq. ( |Q^ ) 
gives the desired boundary condition (|A.18|) . Thus, 5F = if and only if a configuration 
{£, M} satisfies Eq. (|A.17|) with ( |A.18| ), as was to be shown. 

Remark 1: A solution to Eq. (| A . 1 7|) is dynamically stable if the free energy assumes a 
minimum. 

Remark 2: Our Theorem 1 is a finite-temperature generalization of the result obtained for 



cold, catalyzed matter |26 



17 



References 

[1] J.R. Oppenheimer and G.M. Volkoff, Phys. Rev. 55 (1939) 374. 
[2] R. D. Viollier, Prog. Part. Nucl. Phys. 32 (1994) 51. 

[3] J. L. Tonry, Astrophys. J. Lett. 283 (1984) 37; Astrophys. J. 322 (1987) 632. 
[4] A. Dressier, Astrophys. J. 286 (1994) 97. 

[5] A. Dressier and D. O. Richstone, Astrophys. J. 324 (1988) 701. 

[6] J. Kormendy, Astrophys. J. 325 (1988) 128. 

[7] J. Kormendy Astrophys. J. 335 (1988) 40. 

[8] J. Kormendy and D. Richstone, Astrophys. J. 393 (1992) 559. 

[9] J. H. Lacy, C. H. Townes and D. J. Hollenbach, Astrophys. J. 262 (1982) 120. 
[10] J. H. Lacy, J. M. Achtermann and E. Serabyn, Astrophys. J. Lett. 380 (1991) 71. 
[11] D. Lynden-Bell, Nature 223 (1969) 690. 

[12] Ya. B. Zel'dovich and I. D. Novikov, "Stars and Relativity" in Relativistic Astrophysics, 
Vol.1 (U. of Chicago Press, Chicago, 1971). 

[13] R. D. Blandford and M. J. Rees, Mon. Not. R. Astron. Soc. 169 (1974) 395. 

[14] M. C. Begelman, R. D. Blandford and M. J. Rees, Rev. Mod. Phys. 56 (1984) 255. 

[15] N. Bilic, D. Tsiklauri, and R.D. Viollier, Prog. Part. Nucl. Phys. 40 (1998) 17. 

[16] D. Tsiklauri and R.D. Viollier, Astrophys. J. 500 (1998) 591. 

[17] W. Thirring, Z. Phys. 235 (1970) 339. 

[18] P. Hertel and W. Thirring, Comm. Math. Phys. 24 (1971) 22; P. Hertel and W. Thirring, 
"Thermodynamic Instability of a System of Gravitating Fermions", in Quanten unci 
Felder, edited by H. P. Diirr (Vieweg, Braunschweig, 1971). 

[19] P. Hertel, H. Narnhofer and W. Thirring, Comm. Math. Phys. 28 (1972) 159. 

[20] B. Baumgartner, Comm. Math. Phys. 48 (1976) 207. 

[21] J. Messer, J. Math. Phys. 22 (1981) 2910. 

[22] N. Bilic and R.D. Viollier, Phys. Lett. B 408 (1997) 75; N. Bilic and R.D. Viollier, Nucl. 
Phys. B (Proc. Suppl.) 66 (1998) 256. 



18 



[23] J. Ehlers, "Survey of General Relativity Theory" , in Relativity, Astrophysics and Cos- 
mology, edited by W. Israel (D. Reidel Publishing Company, Dordrecht/Boston, 1973), 
sect. 3. 



[24 
[25 
[26 

[27 
[28 
[29 
[30 
[31 
[32 

[33 
[34; 
[35 
[36 
[37; 



J.M. Bardeen, B. Carter, and S.W. Hawking Commun. Math. Phys. 31 (1973) 161. 

A. D. Rendall and B.G. Schmidt, Class. Quantum Grav. 8 (1991) 985. 

B. K. Harrison, K.S. Thorne, M. Wakano and J. A. Wheeler, Gravitation Theory and 
Gravitational Collapse (The University of Chicago Press, Chicago, 1965). 

R.P Feynman, N. Metropolis, and E. Teller, Phys. Rev. 75 (1949) 1561. 

R.F. Sawyer and D.J. Scalapino, Phys. Rev. D 7 (1972) 953. 

E. Witten, Phys. Rev. D 30 (1984) 272. 

W.Y. Chau, K. Lake, and J. Stone, Astrophys. J. 281 (1984) 560. 

A. Kull, R.A. Treumann, and H. Bohringer, Astrophys. J. Lett. 466 (1996) 1. 

N. Bilic, F. Munyaneza, and R.D. Viollier, "Stars and Halos of Degenerate Heavy- 
Neutrino and Neutralino Matter", |astro-ph/9801262| , submitted to Phys. Rev. D. 

F. H Shu, Astrophys. J. 225 (1978) 83. 

P.-H. Chavanis and J. Sommeria, Mon. Not. R. Astron. Soc. 296 (1998) 569. 

L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon, Oxford, 1959). 

R.C. Tolman, Relativity Thermodynamics and Cosmology (Clarendon, Oxford, 1934). 

G. W. Gibbons and S.W. Hawking, Phys. Rev. D 15 (1977) 2752. 



19 



